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Azimuthal distributions of rj and n mesons emitted at 
midrapidity in collisions of 1.9 AGeV 58 Ni+ 58 Ni and 2 AGeV 
40 Ga+ nat C& are studied as a function of the number of 
projectile-like spectator nucleons. The observed anisotropy 
corresponds to a negative elliptic flow signal for r\ mesons, 
indicating a preferred emission perpendicular to the reaction 
plane. The effect is smallest in peripheral Ni+Ni collisions. 
In contrast, for tt° mesons, elliptic flow is observed only in pe- 
ripheral Ni+Ni collisions, changing from positive to negative 
sign with increasing pion transverse momentum. 

PACS numbers: 25.75.-q; 25.75.Dw; 25.75.Ld 

Significant compression of nuclear matter achieved 
during relativistic heavy-ion collisions manifests itself 
in anisotropies of observed azimuthal distributions of 
baryons emerging from the collisions. In-plane emis- 
sion of baryons in the direction of projectile-like (for- 
ward hemisphere) or target-like (backward hemisphere) 
spectator nucleons is designated as positive directed flow, 
while out-of-plane emission of baryons at midrapidity is 
designated as negative elliptic flow, see jlj. 
Similar to baryons, negative elliptic flow was observed 
for high transverse momentum neutral and charged pi- 
ons emitted at midrapidity in 1 AGeV Au+Au collisions 
at SIS (GSI) pflj- I n contrast to positive directed flow 
of baryons, directed flow of charged pions was found to 
be negative The origin of these anisotropies has 

not been attributed to the expansion of nuclear matter, 
rather to the final state interactions of pions with the 
spectator matter transiently concentrated in the reac- 
tion plane P^]. Therefore, an increase in the magnitude 
of elliptic flow of pions with increasing size of spectator 
matter was predicted ||. However, recent detailed study 
of elliptic flow of charged pions in Bi+Bi collisions at 
0.4-1.0 AGeV reports a very small variation of its mag- 
nitude with the change of the number of spectator nucle- 
ons, see g. Hence, the origin of the observed azimuthal 



anisotropies of pions is still unclear. 
While the main source of pions in relativistic heavy-ion 
collisions is the decay of A(1232) resonances, 77 mesons 
are produced mainly by the decay of N*(1535) reso- 
nances. The absorption of pions in hot nuclear matter 
proceeds via the A(1232) resonance which decays domi- 
nantly by pion emission. However, only about 50% of all 
N*(1535) resonances created by 77-meson absorption will 
rccmit rj mesons. Therefore, a comparison of the 77- and 
7T° azimuthal anisotropy may yield information on the 
propagation of these mesons, as well as on the dynamics 
of the parent baryon resonances. 

Below we present results of the first experimental study 
of azimuthal distributions of 77 mesons emitted in colli- 
sions of 1.9 AGeV 58 Ni+ 58 Ni and 2 AGeV 40 Ca+" Qt Ca 
nuclei, and compare them with azimuthal distributions 
of 7T° mesons in the same colliding systems. The experi- 
ments were performed at the Heavy-Ion Synchrotron SIS 
at GSI Darmstadt. 

In the first experiment a 1.9 AGeV 58 Ni beam with an 
intensity of 6.5 x 10 6 particles per spill (spill duration 8 s 
and repetition rate 15 s) was incident on a 58 Ni target 
(502 mg/cm 2 ). In the second experiment, a ™ at Ca target 
(320 mg/cm 2 ) was bombarded by a 40 Ca beam with ki- 
netic energy 2 AGeV and an intensity of 5 x 10 6 particles 
per spill (spill duration 10 s and repetition rate 14 s). 
Photon pairs from the neutral-meson decay were detected 
in the Two- Arm Photon Spectrometer (TAPS) [§. This 
detector system consisted of 384 BaF2 scintillators ar- 
ranged in 6 blocks of 64 modules with individual Charged 
Particle Veto detectors (CPV) in front of each module. 
The blocks were mounted in two towers positioned at 
40° with respect to the beam direction at the distance of 
150 cm. Three blocks were positioned in each tower at 
+21°, 0° and -21° with respect to the horizontal plane. 
In this setup, only neutral mesons around mid-rapidity 
Vcm were detected. The geometrical acceptance of TAPS 
for the 7T and 77 detection was of an order 1 x 10~ 3 . 
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The reaction centrality was determined by the hit multi- 
plicity of charged particles (M react ) in a reaction detec- 
tor. This detector, comprising 40 small plastic scintilla- 
tors, was positioned close to the target and covered the 
polar angles from 14° to 30°. Most of the particles emit- 
ted in this angular range are participant nucleons . 
The plastic Forward Wall (FW) of the KaoS collabora- 
tion, sec Q , comprising 320 plastic scintillators was po- 
sitioned 520 cm downstream of the target and covered the 
polar angles from 0.7° to 10.5°. Particles emitted in this 
angular range are predominantly projectile-like spectator 
nucleons bounced off in the reaction plane Jl(J] . The FW 
provided the information on emission angle, charge and 
time-of-flight of protons and light charged fragments up 
to Z=8. 

The total charge Zpw of particles detected by the FW 
allowed us to estimate the average number of projectile- 
like spectators (A sp ) for each studied bin in multiplicity 
M react , determined by the reaction detector. We used 
the relation (A sp ) = (Z fw) A pro j / 'Z pro j , where A pro j and 
Z P roj are mass and charge of the projectile, respectively. 
The resulting distributions of the total charge Zpw for 
studied bins in M reQCt are shown in the left column of 
Fig.l, both for the Ni+Ni and Ca+Ca collisions. The 
systematic error of the values (A sp ) was found to be less 
than 4 units, see ]l2]j . 

The reaction plane for each event was defined by the in- 
cident beam direction and the vector 



Q = ZiXij | Xi 



where the sum runs over all n particles detected by the 
FW in the event. Xi is the position vector of particle i 
in the x-y plane in the FW and Zi is its charge. The 
unbiased azimuthal-angle distribution of the vector Q 
was found to be flat after introducing a correction, which 
takes into account the shift of the beam position with re- 
spect to the geometrical center of the FW. The required 
corrections were below 1 cm. 

Because of finite multiplicity fluctuations, the azimuthal 
angle (f> of the vector Q can differ from the azimuthal an- 
gle of the true reaction plane fytme by a deviation A(j) p i. 
To estimate the width a p i of the distribution N(A<j) p i) 
we randomly divided the hits in each event into two sub- 
groups containing each one half of the number of parti- 
cles, see |fL3[ . For each subevent one can construct two 
independent vectors Q\ and Q 2 , according to Eq.(l) and 
extract the angle Aci>i2=0i —<h between the two vectors. 
The width o\ 2 of the resulting distribution N(A<fii2) is 
a measure of the precision of the reaction-plane determi- 
nation. In the approximation of a Gaussian distribution 
of N(A(f>pi) one has a p i = 1712/2. Our analysis was re- 
stricted to a sufficient vector length | Q | in order to re- 
ject the most central events with no spectator flow. This 
selection removes 25% of all registered events. The reso- 



lution a p i varies between 43° and 55° depending on the 
reaction centrality and the colliding system, see Table Q. 
These values agree with published data from studies of 
charged-baryon flow in similar colliding systems [p^|-p^[ . 
Photon-particle discrimination has been performed in 



TAPS as described in 16 1. For each pair of detected pho- 
tons in a given event we calculated the invariant mass 
M pa i r and the momenta p pa i r using the following rela- 
tions: M^ air = 2E 1 E 2 (l - cosQ 12 ) and p pair = pi + p 2 , 
where E\, pi and E 2 , p 2 are the energies and momenta 
of the corresponding photons, 812 is the opening angle 
of the photon pair. We analyzed only neutral mesons 
in a narrow rapidity window (\y — y cm \ < 0.1). The 
combinatorial background, due to uncorrelated photon 
pairs, was deduced by the method of event mixing from 
experimental data, i.e. combining photons from different 
events with the same overall event characteristics ||Tci[ , [L6| . 
The resolution in M pair was 15 MeV/c 2 and 45 MeV/c 2 
(FWHM) for the ir° and 77 peaks, respectively fl5| . 
We found that only the magnitude of the combinato- 
rial background is dependent on the azimuthal angle 
Aip = (ppair — (j) of meson emission relative to the reac- 
tion plane Therefore, the shape of the combinatorial 
background deduced by the method of event mixing for 
each bin in M react and transverse momenta pt of the pho- 
ton pair was scaled to the experimental data for each bin 
in azimuthal angle Aip separately and subtracted from 
the data. We present in Fig.l the resulting azimuthal 
yields of r\ mesons, soft ir° ( 0< pt < 200 MeV/c ) and 
(X) hard tt° ( 600< p t < 800 MeV/c ) for different bins in 
M react for both systems studied. 

We fitted the azimuthal yields of rj and ir° mesons both by 
assuming a constant ( isotropic distribution ) and by the 
first two terms of a Fourier expansion in the azimuthal 
angle: 



iVn 

N(Atp) = — - ( 1 + 2v lC os(Aip) + 2v 2 cos(2A<p) ) , (2) 
2n 

which are used to parametrize the directed (vi) and el- 
liptic (v 2 ) flow [0. The standard F-test rejects the fit 
of azimuthal yields of 77 mesons by a constant with a con- 
fidence level of 95%. 

The extracted values of v\ are zero within the error bars, 
as should be expected since we study symmetric colliding 
systems at midrapidity, see EH]. 

The resulting values of the parameter v 2 for 77 and 7T° 
mesons are given in Table Q. The parameter v 2 is nega- 
tive for -q mesons, indicating a preferred emission perpen- 
dicular to the reaction plane (negative elliptic flow). In 
contrast to previous studies of heavy systems at 1 AGeV 
the elliptic-flow signal v 2 for tt° mesons is pronounced in 
peripheral {(A sp ) = 51) Ni+Ni collisions only. 
The measured azimuthal distributions are affected by the 
resolution in the determination of the reaction plane. 
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FIG. 1. The left part shows the measured distribution of the total charge of particles detected in the FW for different bins 
in multiplicity M react (see Table |): (a) - (c) for the experiment 58 Ni+ 58 Ni at 1.9 AGeV and (d) - (e) for the experiment 
40 Ca+" at Ca at 2 AGeV. The arrows indicate mean values of total charge. The right part shows the measured azimuthal angle 
distributions of r\ mesons, soft ty° (0< pt < 200 MeV/c) and hard 7r° (600< pt < 800 MeV/c) with respect to the reaction plane 
which correspond to the same bins in reaction centrality (see left part of the picture). The horizontal error bars indicate the 
bin width. 



Therefore, it is necessary to correct the V2 coefficients 
for the fluctuation Acf> p i = (j>true — 4> oi the azimuthal 
angle of the estimated reaction plane with respect to the 
true one. Averaging over many events, one obtains the 
following relation between the measured V2 coefficients 
and the true w^" 6 coefficients: v l 2 U(L — v 2 / (cos2A<f) p i) 
[p~3| , P~r|| . In our analysis we used the method described in 
detail in Q, which allows to directly extract the values 
(cos2A(fipi) by means of an analytical expression from the 
experimental distribution of A r (A0i2) (see above). The 
obtained values of {cos2A<j) p i) are given in Table | for 
each bin in reaction centrality. The systematic error of 
the values (cos2A0 p /) was found to be less than 15% [p"5[ . 
The dependence of v% ue on the average number of 
projectile-like spectator nucleons (A sp ) is shown in Fig. 2 
for rj mesons and for two bins in p t of tt° mesons for both 
Ni+Ni and Ca+Ca collisions. 

As it was already mentioned above, microscopic calcula- 
tions associate elliptic flow of pions with the assumption 
of their "shadowing" by spectators. Guided by this we 
attempt to describe our data by a geometrical model as- 
suming absorption of the mesons in the spectator blobs. 



We calculate the parameter v| r " e as function of the num- 
ber of spectators A sp from the equation 



= 0.5(1 -R)/(l + R) 



(3) 



with R — exp(L/X), where A is the mean free path for 

1 /3 

mesons in cold spectator matter and L = 2 • A sp fm is 
the mean thickness of spectator matter. The mean free 
path in cold spectator matter for rj mesons is A n ~ 1-2 fm 
for the momentum range p v w 50-200 MeV/c fl^ , |l9|| , and 
for 7T° mesons A ff o w 1-6 fm |p0|| . The dash-dotted lines 
in Fig. 2 present the results of these calculations for two 
different values of A: X—2 fm and A=6 fm. It is obvious, 
that the simple scenario assuming final-state interactions 
with spectators alone can not describe the strong differ- 
ence between the rj and tt° azimuthal anisotropics. The 
observed dependence of the magnitude of the azimuthal 
anisotropy for rj mesons on the number of spectator nu- 
cleons seems to contradict the model, too. 
The rj mesons are solely produced by the decay of the 
heavy N*(1535) resonances which can only be excited in 
the early stage of the collision. Consequently, most of the 
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r\ mesons are emitted after rescattering in the early phase 
of the collision, while most of the 7r° mesons are emitted 
after the spectators have passed the collision zone. This 
may explain the much stronger azimuthal anisotropy ob- 
served for 77 mesons. However, for a quantitative expla- 
nation microscopic model calculations are needed. 
« 0.4 1 1 
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FIG. 2. The parameters v^ ue for rj mesons (squares), 
soft 7T° ( 0< p t < 200 MeV/c ) (triangles) and hard n° 
( 600< p t < 800 MeV/c ) (circles) as a function of the 
number of spectators (A sp ). Full symbols correspond to 
Ca+Ca and open symbols correspond to Ni+Ni collisions. 
The dash-dotted lines represent results of the geometrical ab- 
sorption model (see text). 

In summary, we have studied simultaneously azimuthal 
angular distributions of ir° mesons and of 77 mesons emit- 
ted at midrapidity in the two colliding systems 58 Ni+ 58 Ni 
at 1.9 AGeV and 40 Ca+" at Ca at 2 AGeV. We observed 
a strong out-of-plane elliptic flow of rj mesons. The el- 
liptic flow of 7T° mesons is very weak in contrast to data 
obtained for heavy colliding systems at 1 AGeV. 
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TABLE I. Parameters V2 (not corrected for the reaction-plane uncertainty) for ir° and rj mesons deduced from the experi- 
mental azimuthal distributions for several intervals (a-e) in M reac t. Corresponding mean number of projectile-like spectators 
(Asp), uncertainty of the reaction plane angle a p i, correction of V2 due to this uncertainty are given (see text). For 7r° mesons 
the selected intervals in transverse momentum p t are indicated. 
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58 Nii 58 Ni at 1.9 AGeV 




40 Cai nat 


Ca at 2 AGeV 




^1 reac 


(a) 2 - 6 


(b) 7 - 10 


(c) > 11 


(d) 1 - 3 


(e) > 4 




(A sp ) 


51 


44 


37 


34 


28 




ffpi [°] 


43 


46 


50 


52 


55 
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cos2A(p p i) 


0.37 


0.32 


0.26 


0.22 


0.19 




pt MeV/c 












v 2 (rj) 


- 600 


-0.08T0.04 


-0.12i0.05 


-0.10i0.06 


-0.09i0.03 


-0.09i0.04 




) 0-200 


0.033T0.017 


0.021i0.019 


0.026i0.031 


0.012i0.016 


-0.006i0.021 




200-400 


0.005±0.007 


0.007i0.009 


0.011i0.016 


-0.011i0.006 


-0.005i0.011 




400-600 


-0.004±0.005 


0.008i0.011 


-0.018i0.016 


0.006i0.007 


0.016i0.007 




600-800 


-0.031T0.012 


0.012i0.013 


0.022i0.021 


0.005i0.012 


0.007i0.012 



4 



